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1  Summary 


The  research  supported  under  this  grant  focussed  on  the  accurate  com¬ 
putation  of  compressible  leading-edge  vortex  flows  with  interpretation  and 
understanding  of  the  results.  To  this  end,  a  number  of  investigations  were 
undertaken,  starting  with  the  conical  Euler  equation  model  and  concluding 
with  fully  three-dimensional  Navier-Stokes  equation  calculations.  Consider¬ 
able  emphasis  was  placed  on  . 

♦  Comparisons  of  Euler  with  Navier-Stokes  calculations  to  understand 
the  origins  of  various  mechanisms^ 

r*  Comparison  of  calculations  with  experiments  both  to  determine  the 
accuracy  of  the  calculations  and  to  comprehend  the  physical  processes’ 
and 

V  Interpretation  of  the  results  with  the  aid  of  analytical  models. 

A  summary  of  some  of  the  work  is  contained  in  the  paper  by  Murman, 
Goodsell  ,  and  Malecki  [1]  presented  at  Symposium  Transsonicum  III  which 
is  attached  as  Appendix  I.  The  major  contributions  are  contained  in  the 
thesis  of  the  participating  students,  namely: 

•  Doctoral  Students 


1.  Kennth  Grant  Powell,  Vortical-Solutions  of  the  Conical  Euler 
Equations  [2] 

2.  Bernard  Loyd,  A  Semi-Implicit  Navier-Stokes  Solver  and  its  ap¬ 
plication  to  a  Study  of  Separated  Flow  about  Blunt  Delta  Wings 

[3] 

•  Masters  Students 

1.  Robert  E.  Malecki,  Euler-Equation  Calculations  for  a  Cropped- 
Delta  Wing  Using  the  CYBER  205  [4] 

2.  Aga  Myung  Goodsell,  3-D  Euler  Calculations  of  Vortex  Flows 
Over  Delta  Wings  [5] 

3.  Jorge  P6rez-S4nchez,  Numerical  Simulation  of  Deceleration  of  an 
Axisymmetric  Vortex  [6]. 


Other  papers  and  articles  resulting  from  the  research  are  cited  in  the  fol¬ 
lowing  sections  and  listed  in  the  References  section. 


2  Conical  Euler  Solutions 


The  conical  Euler  equation  model  was  adopted  to  provide  a  computation¬ 
ally  tractable  approach  for  understanding  the  essential  features  of  vortical 
flow  represented  by  inviscid  mechanisms.  Conical  calculations  are  valid  for 
steady  supersonic  freestream  Mach  numbers  and  bodies  which  have  conical 
self-similarity.  The  wings  must  have  a  sharp  leading  edge  so  that  the  Kutta 
condition  yields  a  realistic  model  for  the  separation  point. 

The  PhD  thesis  of  Powell  [2]  contains  a  thorough  description  of  the  nu¬ 
merical  method,  numerous  computed  results,  comparisons  with  experiments, 
and  interpretation  of  the  findings.  Calculations  were  done  for  delta  wings 
at  various  angles  of  attack  and  yaw,  and  for  delta  wings  with  leading  edge 
flaps.  It  was  found  that  to  achieve  adequate  mesh  resolution  in  the  region  of 
the  leading  edge  vortex,  local  refinement  was  needed.  An  embedded  mesh 
strategy  was  developed  and  applied  to  most  of  these  cases.  This  is  reported 
in  the  paper  by  Powell  and  Murman  [7].  Figures  1-4  on  the  following  pages 
present  representative  results  from  this  study.  These  are  for  a  flat  plate 
delta  wing  with  a  75°  leading  edge  sweep  at  12°  angle  of  attack,  8°  of  yaw 
and  a  freestream  Mach  number  of  1.7. 

A  parallel  investigation  [8]  involved  a  comparsion  of  the  conical  Euler 
equaton  results  with  experiments  and  conical  Navier-Stokes  equation  cal¬ 
culations.  The  cases  were  all  for  supersonic  freestream  Mach  numbers,  but 
included  three  leading  edge  sweep  angles  and  five  angles  of  attack.  Although 
the  Euler  equation  solutions  predicted  the  major  flow  regimes  quite  well,  the 
Navier-Stokes  solutions  were  in  somewhat  better  agreement  with  the  data 
as  might  be  expected.  This  paper  was  selected  as  the  best  Applied  Aerody¬ 
namics  Paper  in  1987  by  the  ALAA. 
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Figure  1  -  Grid  showing  two  levels  of  embedding 
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2  -  Cross-flow  streamlines  showing  topology  of  the  vortices 
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Figure  4  -  Comparison  of  surface  pressure  coefficient  with  data 
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3  Total  Pressure  Loss 


A  question  of  considerable  interest  is  what  is  the  origin  of  the  large  total 
pressure  losses  in  the  cores  of  the  vortices,  as  indicated  in  Figure  3.  The 
Euler  equations  alone  do  not  contain  any  dissipative  mechanism.  However, 
in  their  discrete  form,  both  added  and  implicit  artifical  viscosity  are  present 
and  necessary.  A  paper  was  written  by  Murman,  et  al  [9]  which  discussed 
this  issue.  The  paper  contained  an  error  in  the  cross-flow  streamline  integra¬ 
tion  scheme  which  led  to  some  wrong  conclusions  that  were  later  corrected 
[!°]. 

The  work  of  Powell  [2]  conclusively  showed  that  the  total  pressure  losses 
are  by  and  large  independent  of  both  second  and  fourth  difference  damping 
coefficients,  (j  and  £4  respecitvely,  as  well  as  grid  refinement.  The  tables 
on  the  following  page  present  these  results.  Figures  5  and  6  and  the  tables 
indicate  that  the  overall  vortex  dimensions  (width  and  height)  are  less  sen¬ 
sitive  to  these  parameters  than  is  the  region  of  large  total  pressure  loss  (core 
size). 

Powell  proceeded  to  formulate  a  theoretical  model  for  the  core  of  a  vis¬ 
cous  vortex  by  deriving  a  self-similarity  theory  for  an  incompressible,  conical 
vortex  [2,11].  The  model  predicted  that  the  total  pressure  loss  is  determined 
by  the  magnitude  of  the  swirl  velocity  at  the  edge  of  the  vortex,  and  not  by 
the  magnitude  of  the  viscosity.  The  total  pressure  loss  is  related  to  the  dis¬ 
sipation  mechanism,  and  apparently  the  numerical  viscosity  is  a  reasonable 
representation  of  physical  viscosity  for  some  cases.  Powell  also  showed  that 
the  Reynolds  number  based  upon  artifical  viscosity  is  so  large  in  the  region 
of  the  vortex  that  it  would  dominate  real  visous  terms  in  most  Navier-Stokes 
calculations.  This  is  discussed  further  in  Appendix  I. 

Finally,  based  upon  the  observation  that  the  total  pressure  loss  in  the 
vortex  might  be  related  to  bursting,  a  preliminary  study  was  done  to  look 
at  the  effects  of  decelerating  an  axisymmetric  inviscid  rotational  flow  [6]. 
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Loss  Level 

Core  Size 

Vortex  Width 

Vortex  Height 

0.003 

0.0010 

37.4% 

1.00 

1.00 

1.00 

0.003 

0.0003 

37.6% 

0.32 

1.10 

1.13 

0.003 

0.0030 

38.2% 

1.74 

1.08 

1.02 

0.001 

0.0010 

37.2% 

0.71 

0.97 

1.01 

0.010 

0.0010 

37.6% 

1.67 

1.11 

1.08 

Table  Effects  of  artificial  viscosity  level  on  loss 


Equivalent  grid 

Loss  Level 

Core  Size 

Vortex  Width 

Vortex  Height 

64  x  64 

33.9%' 

2.55 

1.39 

1.43 

128  x  128 

37.6% 

1.00 

1.00 

1.00 

256  x  256 

35.1% 

0.26 

0.96 

1.05 

Table  2:  Effects  of  grid  refinement  on  loss 
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Figures  5  and  6  -  Total  pressure  losses;  64  x  64  Sc  256  x  256  grids. 
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4  Three-Dimensional  Euler  Equations 

Three  studies  were  undertaken  [4,5,12]  based  upon  three-dimensional 
Euler  equation  calculations.  In  the  first  study,  Malecki  [4]  vectorized  an 
existing  Euler  solver  for  the  CYBER  205  and  used  an  0-0  grid  supplied  by 
Rizzi  to  compute  the  flow  about  a  cropped  delta  wing  model  tested  in  the 
Netherlands.  The  study  showed  that  significant  levels  of  artificial  viscosity 
were  needed  to  obtain  an  converged  solution,  apparently  due  to  high  grid 
skewness. 

Goodsell  [5]  revamped  the  solver  to  perform  on  an  O-H  grid  and  pro¬ 
ceeded  to  compute  not  only  the  same  transonic  cases  that  Malecki  had  done, 
but  also  supersonic  and  subsonic  cases.  Representative  results  of  both  of 
these  studies  are  contained  in  the  paper  of  Appendix  I  [1], 

The  third  study  done  by  Lee  and  Murman  [12]  utilized  Goodsell’s  Euler 
solver  for  a  preliminary  study  of  the  effects  of  the  vortex  from  a  closely 
coupled  canard  on  the  following  wing.  Representative  results  are  given  in 
Figures  7  and  8. 

The  three-dimensional  Euler  calculations  showed  that  the  absence  of  sec¬ 
ondary  separation  which  is  not  included  in  the  Euler  model  had  substantial 
effects  on  the  predicted  pressure  distributions  in  the  transonic  regime,  ne¬ 
cessitating  the  development  of  a  Navier-Stokes  model.  Supersonic  cases  are 
much  less  sensitive  to  this  effect  [8]. 
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Figure  7:  Canard-on  pressure  surface  distribution 
a.  Computational  b.  Experimental 
a  =  8.8* 

® 


Y  as  1.2 
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Figure  8:  Canard-off  surface  pressure  distribution 
a.  Computational,  b.  Experimental 
a  =  8.8* 


5  Three-Dimensional  Navier-Stokes  Solutions 


Two  studies  were  done  for  three-dimensional  Navier-Stokes  equations. 
The  first  was  undertaken  by  Rizzi  who  spent  a  summer  at  MIT  as  part  of 
this  project.  In  this  regard,  he  modified  his  existing  Euler  solver  to  add 
the  full  Navier-Stokes  terms,  and  then  proceeded  to  compare  the  Euler  and 
Navier-Stokes  results  [13]  for  a  round  leading  edge  version  of  the  cropped 
delta  wing  studied  by  Malecki  and  Goodsell.  Representative  results  shown  in 
Figure  9  indicate  that  the  viscous  calculations  are  in  much  better  agreement 
with  experiments  than  are  the  inviscid  ones. 

The  second  study  was  undertaken  by  Loyd  [3,14]  who  added  the  thin 
layer  Navier-Stokes  terms  to  Goodsell’s  code  and  used  a  newly  developed 
semi-implicit  method  to  solve  them.  Loyd  particularly  focussed  on  the  mech¬ 
anisms  of  leading  edge  separation  for  laminar  flows  over  a  family  of  swept 
conical  delta  wings.  Detailed  analysis  indicated  both  shock-free  and  shock- 
induced  separation  at  the  leading  edge.  The  other  features  of  primary  and 
secondary  vortical  flows  were  also  observed.  A  representative  result  from 
this  study  is  shown  in  Figure  10. 


9 


x/cs0.6 


x/cs0.8 


Figure  9  -  Measured  and  computed  C-  distributions  [13] 


Figure  10  -  Numerical  oil  flow  picture,  =  1.6,  a  =  4' 
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Transsonic  Vortical  Flow 
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Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


Introduction 

Transonic  flow  is  generally  associated  with  the  inviscid  fluid  effects  when  a  flow  with  a 
freestream  Mach  number  M^  ~  1  accelerates  to  locally  supersonic  velocities,  or  decelerates 
to  locally  subsonic  velocities,  as  it  moves  streamwise  past  a  body.  However,  transonic  flows 
can  occur  in  many  other  circumstances  such  as  on  high  lift  devices  at  low  Moo,  at  the  tips  of 
helicopter  blades,  and  around  corners  in  ducts.  In  this  paper,  a  class  of  transonic  flows  for 
slender,  highly  swept  wings  at  high  angle  of  attack  (a  >  10°)  will  be  considered.  Under  these 
conditions,  leading  edge  vortices  generally  form  as  shown  schematically  in  Figure  1.  If  Moo  ~  1, 
the  streamwise  flow  can  be  transonic  in  the  classical  sense,  leading  to  mixed  subsonic-supersonic 
flow  and  shock  waves  on  the  wing.  However,  the  leading  edge  vortices  create  additional  transonic 
cross  flow  effects  due  to  the  high  swirl  velocities  associated  with  them. 


Figure  1:  Schematic  Flow  Field  Figure  2:  Classification  of  Flow  Regimes 

The  classes  of  vortical  transonic  crossflows  are  generally  correlated  in  terms  of  the  angle  of 
attack  and  Mach  number  normal  to  the  wing  leading  edge  given  by  the  formulas  [1] 

On  —  tan-1  (tan  a  sec  A)  (1) 

Mn  =  Moo  sin  a  esc  On  (2) 

where  A  is  the  wing  leading  edge  sweep  angle.  Figure  2  shows  a  correlation  from  experimental 
data  presented  and  discussed  by  Miller  and  Wood  [2,3]  for  sharp  leading  edge  wings.  Other  flow 
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’Present  Address,  Pratt  it  Whitney  Aircraft  Group,  Hartford,  Conn. 


patterns  have  also  been  discovered  experimentally  by  Vorropoulos  and  Wendt  [4]  and  Szodruch 
and  Peake  {5j.  A  paper  in  this  conference  by  Narayan  and  Seshadri  [6]  further  discusses  these 
Sow  regimes  and  also  presents  results  for  round  leading  edge  wings.  For  the  latter  case,  a 
Reynolds  number  based  on  some  leading  edge  length  scale  must  also  be  an  important  parameter. 

Transonic  vortical  Sows  are  a  complex  mixture  of  streamwise  and  crossflow  effects.  The 
paper  by  Elsenaar  and  BiiteSsch  [7]  presents  results  of  an  extensive  experimental  investigation 
for  a  65s  delta  wing  which  serve  to  illustrate  this  point.  A  complete  model  which  can  handle 
the  interacting  inviscid  and  viscous  effects  undoubtedly  requires  Navier-Stokes  equations,  and 
one  such  approach  is  given  by  Rizzi,  Drougge,  and  Muller  [8].  For  sharp  leading  edge  wings,  the 
separation  of  the  flow  to  form  a  primary  vortex  is  fixed  by  the  geometry  via  a  Kutta  condition, 
and  an  Euler  equation  model  should  be  adequate  if  vorticity  is  introduced.  However,  important 
viscous  effects  such  as  shock  wave  boundary  layer  interaction,  secondary  vortex  formation,  and 
boundary  layer  transition  clearly  won’t  be  represented.  Although  the  question  of  whether  such 
models  can  represent  vortex  breakdown  is  not  settled,  numerical  solutions  seem  to  indicate 
that  bursting  of  the  primary  vortex  is  captured  if  the  mesh  is  fine  enough.  In  principle,  a 
potential  flow  with  added  vortex  sheets  or  filaments  as  suggested  by  Desai  et  ad  [9]  could  be 
used.  However,  the  direct  solution  of  the  Euler  and/or  Navier-Stokes  equations  seems  more 
promising,  although  more  expensive. 

In  this  paper,  Euler  equation  models  are  considered  for  sharp  edge  geometries.  The  conicad 
Euler  equations  are  considered  first  to  investigate  the  crossflow  phenomena.  Following  this 
is  a  discussion  of  the  origin  of  the  vorticity  in  Euler  equation  calculations.  To  conclude  the 
paper,  three-dimensional  solutions  are  presented  which  combine  both  streamwise  and  crossflow 
transonic  effects.  The  geometry  considered  is  that  reported  by  Elsenaar  and  Biitefisch  [7]. 

Governing  Equations 

The  three-dimensional  unsteady  Euler  equations  in  conservation  form  are: 
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(u,  v,  w)  =  u  are  the  Cartesian  velocity  components  in  the  z,  y,  z  directions,  respectively,  p  is 
the  density,  E  the  total  energy  per  unit  mass,  and  p  the  static  pressure  given  by  the  equation 
of  state 
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Unsteady  terms  have  been  retained  for  the  iterative  procedure  used  to  reach  the  steady  solution. 
The  tangential  flow  boundary  condition  is  enforced  on  the  solid  body  while  the  Kutta  condition 
will  be  enforced  implicitly  by  the  numerical  method.  Freestream  conditions  are  applied  far 
away  from  the  body. 

Introducing  the  conical  variables 


r  =  Jz2  +  yJ  +  z 1 


and  assuming  conical  self-similarity  (the  solution  is  independent  of  r),  the  Euler  equations 
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These  equations  are  valid  for  wing  geometries  generated  by  rays  passing  through  the  origin  and 
for  A/*  >  1  so  that  the  outer  boundary  condition  is  also  conical.  Again,  the  unsteady  terms  are 
retained  only  for  the  iterative  solution  process.  The  equations  are  solved  on  the  unit  sphere  by 
setting  r  =  1.  Boundary  conditions  are  the  same  as  the  Euler  equations  except  that  freestream 
conditions  are  applied  ahead  of  the  conical  bow  shock. 


Solution  Procedure 


The  basic  solution  scheme  is  a  finite  volume  spatial  discretization  with  a  multi-stage  inte¬ 
gration  in  time,  as  proposed  by  Jameson  et  a 1  [10].  It  will  be  outlined  for  the  conical  Euler 
equations.  The  solution  of  the  three-dimensional  Euler  equations,  which  is  very  similar,  is  fully 
explained  by  Roberts  [11]  who  implemented  the  version  used  in  this  paper. 


The  conical  Euler  equations  (6)  written  in  vector  form 
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may  be  transformed  using  the  divergence  theorem  into  an  integral  equation  and  discretized  as 


dU, 
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(8) 


Aij  is  the  area  of  a  computational  cell  and  the  summation  represents  the  fluxes  through  the 
four  cell  sides  with  normal  vectors  n*  =  (n^,  nnt).  The  flux  vectors  F*,  G*  on  the  cell  sides  are 
evaluated  from  simple  averages  of  flux  vectors  at  the  cell  centers  F GtJ  and  F,±i,j±i,  Gi±i,j ±i 
.  A  multistage  integration  scheme  is  used  to  solve  the  equations  in  time  to  reach  a  steady 
state.  In  order  to  accelerate  the  iterative  solution,  the  time  step  in  each  cell  is  determined  by 
the  linear  stability  criteria  for  that  cell.  Since  the  time  steps  are  different  for  each  cell,  the 
solution  is  not  time  accurate.  The  tangential  flow  boundary  condition  on  the  body  is  enforced 
by  zeroing  the  convective  fluxes  in  (6)  and  extrapolating  the  pressure  from  the  adjacent  cell 
center,  piu4f  =  pjiml  . 

This  scheme  requires  added  second  and  fourth  order  damping  to  capture  shocks  and  damp 
unwanted  high  frequency  modes  of  the  discrete  solutions.  The  damping  formulation  follows 
Rizzi  and  Eriksson  [12].  The  second  order  damping  is  pressure-weighted,  and  is  of  the  form 


where  and  are  undivided  central  difference  operators  defined  by 

*(U,,)  =  uI+i,  -  UM,  «,( UU)  =  U,,+i  -  U,M. 

The  fourth-order  damping  is  unweighted,  and  of  the  form 

D4(U)  =  vA  [tf  U  +  6/u] . 


(9) 

(10) 

(11) 


The  pressure  switch  in  the  second  order  damping  is  normalized  to  vary  from  O  (A2)  to  1 
throughout  the  flow  field,  where  A  is  a  measure  of  the  mesh  spacing.  It  is  designed  to  turn  on 
at  embedded  shock  waves.  In  these  calculations  it  was  found  that  the  second  difference  damping 
needed  to  be  turned  on  at  the  wing  leading  edge.  This  was  done  by  setting  the  pressure  switch 
to  1  in  a  few  cells  in  this  region.  Values  of  i/j  =  .05  and  t/4  =  .01  were  used  for  the  calculations 
unless  noted. 

The  added  damping  terms  (9)  and  (11)  together  with  the  basic  discretization  error  from  (8) 
implies  that  the  numerical  calculations  represent  the  solution  of  a  modified  differential  equation 

dUJ  d*V' 
d£4  +  dij4 

(12) 

where  5(,5n  are  the  pressure  switches  in  (9)  and  A  the  maximum  characteristic  speed  of  propa¬ 
gation  of  a  disturbance  cs  |uj  +c.  The  first  group  of  terms  on  the  right  hand  side  represents  the 
second  order  damping,  the  second  group  the  truncation  errors,  and  the  third  group  the  fourth 
order  damping.  In  the  neighborhood  of  the  tip  and  embedded  shock  waves,  S(,Sn  ~  1  and  the 
second  difference  damping  is  proportional  to  the  first  power  of  the  mesh  spacing.  Away  from 
these  regions,  ~  A2  and  the  second  difference  damping  is  proportional  to  the  cube  of 

the  mesh  spacing.  The  fourth  difference  damping  is  everywhere  proportional  to  the  cube  of  the 
mesh  spacing.  The  solutions  automatically  satisfy  the  Kutta  condition,  and  the  present  authors 
believe  this  is  due  to  the  added  numerical  viscosity. 


h* + §iG + H  -  *',u  {§(  + I?  (s,^r)] + 0 


Conical  Euler  Solutions 

Numerous  solutions  have  been  computed  for  the  conical  Euler  equations  for  sharp  edge 
geometries  using  the  method  as  outlined  above  or  in  a  somewhat  modified  form.  All  the 
primary  vortex  and  shock  wave  structures  illustrated  in  Figure  2  have  been  reproduced,  as  well 
as  other  shock  structures  reported  in  [4,5].  Secondary  vortices  due  to  boundary  layer  separation 
under  the  primary  vortex  and  vortices  created  by  shock  wave-boundary  layer  interactions  are 
not  modeled  by  the  Euler  equations.  The  most  comprehensive  presentation  of  the  variety  of 
solutions  which  have  been  obtained  may  be  found  in  the  PhD  thesis  of  Powell  [13] .  A  comparison 
of  the  conical  and  a  three-dimensional  Euler  solution  are  given  in  [14],  while  comparisons 
between  the  conical  Euler  and  conical  Navier-Stokes  solutions  with  data  are  given  in  [15].  In 
this  section,  two  representative  examples  will  be  presented. 

An  interesting  example  illustrating  the  nature  of  the  transonic  crossflow  effects  is  a  flat 
plate  delta  wing  with  a  leading  edge  sweep  A  =  75°  in  a  =  1.95  flow  at  an  angle  of  attack 
a  =  25°.  The  above  conditions  correspond  to  Mn  =  94 ,an  =  61°  and  fall  in  the  middle 
top  region  of  Figure  2.  Experimental  results  for  this  condition  are  reported  by  Monnerie  and 
Werle  [16].  The  calculations  were  done  on  an  O-type  grid  which  had  128  cells  around  the  wing, 
windward  to  leeward  symmetry  plane,  and  128  cells  between  the  wing  and  the  farfield  boundary 
outside  the  bow  shock. 

The  general  topology  of  the  flow  field  is  illustrated  in  Figure  3  by  the  conical  streamlines 
which  are  everywhere  tangential  to  the  crossflow  velocity  vectors.  The  view  corresponds  to 
looking  upstream  from  the  trailing  edge  with  only  the  right  hand  side  of  the  wing  shown. 
The  dominant  feature  is  the  sheet  emanating  from  the  wing  leading  edge  and  spiraling  into 
the  vortex  core3.  On  the  windward  side,  the  attachment  line  is  almost  halfway  between  the 
centerline  and  leading  edge  of  the  wing,  while  on  the  leeward  side  it  is  at  the  wing  centerline. 

*The  blank  area  ia  th«  middle  of  the  core  ia  do*  to  aa  inaccurate  itreamllne  integration  [17]  and  aot  aa 
apparent  eource  as  suggeeted  ia  [18]. 


Figure  3:  Conical  streamlines. 


Figure  4:  Total  Pressure  Loss. 


Contours  of  constant  total  pressure  loss  1  —  ( Pt / Pt*,)  &re  shown  in  Figure  4.  There  are  large 
losses  in  the  shear  layer  from  the  leading  edge,  and  in  the  center  of  the  vortex  the  total  pressure 
loss  is  88  %.  It  can  be  seen  that  cross  flow  shocks  are  present  above  and  below  the  primary 
vortex  and  that  there  is  also  a  shock  between  the  two  primary  vortices.  This  latter  shock  is 
observed  in  the  shadowgraph  picture  from  the  experiment  [16].  The  swirling  flow  accelerates  to 
cross  flow  Mach  number  Me  =  2.5  above  the  vortex  before  passing  through  the  first  shock.  It 
then  turns  downwards  towards  the  wing  reaching  Afe  —  1.4  ahead  of  a  second  shock.  The  flow 
then  accelerates  to  Afe  =  1.8  on  the  wing  surface  under  the  vortex  where  it  experiences  the  third 
shock.  A  comparison  of  the  cross  flow  streamlines  and  the  shock  locations  will  reveal  a  sort  of 
transonic  nozzle  effect  for  the  last  two  shocks.  There  is  a  minimum  spacing  between  the  vortex 
and  the  symmetry  plane  and  the  vortex  and  the  wing.  Examination  of  the  crossflow  Mach 
number  contours  which  aren’t  shown  here  reveals  a  converging-diverging  nozzlelike  pattern. 


Figure  5:  Experimental  pitot  pressure.  Figure  6:  Computed  pitot  pressures. 

Monnerie  and  Werl4  reported  measured  pitot  pressure  (normalized  by  the  freestream  total 
pressure)  values  on  the  leeside  of  the  wing  which  are  shown  in  Figure  5.  Pitot  pressures 
computed  from  the  Euler  calculations  using  the  Rayleigh  pitot  formula  (19)  are  shown  in  Figure 
6.  The  computed  values  would  correspond  to  a  probe  which  is  everywhere  aligned  with  the  flow, 
while  in  the  experiments  the  probe  was  aligned  with  the  wing  centerline.  Vortical  flows  have 
large  flow  angularity  as  illustrated  by  the  swirl  angle  (ratio  of  croes  flow  velocity  magnitude  to 


radial  velocity)  contours  of  Figure  7.  Generally,  pitot  pressure  measurements  are  accurate  to 
within  1 %  if  the  angularity  is  less  than  10s  [20] .  This  is  really  only  true  right  in  the  vortex 
core  region,  and  the  computed  (.05)  and  measured  (.06)  values  are  in  excellent  agreement. 
Throughout  the  rest  of  the  region  there  is  good  qualitative  agreement  except  where  the  wing 
boundary  layer  and  secondary  vortex  formation  dominate. 


Figure  7:  Computed  swirl  angles.  Figure  8:  Symmetry  plane  flow  angles. 


Measurements  were  also  made  of  the  flow  angle  in  the  leeside  symmetry  plane.  A  small  wedge 
in  the  flow  was  adjusted  until  the  top  and  bottom  shock  waves  viewed  with  a  shadowgraph  were 
symmetric.  The  estimated  accuracy4  is  ±2°.  Figure  8  compares  the  measured  and  computed 
values  which  are  in  good  agreement.  The  change  in  flow  angularity  across  the  shock  spanning 
the  port  and  starboard  vortices  is  clear. 

Calculations  have  been  done  on  coarser  grids  to  investigate  the  sensitivity  of  the  solution 
to  grid  resolution.  The  following  table  displays  some  of  these  values. 


Grid 

32  x  32 

64  x  64 

128  x  128 

r)  vortex 

.287 

.307 

.320 

£  vortex 

.445 

.470 

.438 

Normal  Force  Coefficient 

.729 

.739 

.742 

Max  Swirl  Angle 

38° 

o 

O 

■w 

e 

O 

■w 

Max  Total  Pretture  lost 

.872 

.882 

.880 

Min  Pitot  Pressure 

.075 

.052 

.048 

It  can  be  seen  that  the  position  of  the  vortex,  the  normal  force  coefficient,  the  maximum  total 
pressure  loss,  and  maximum  swirl  angle  are  all  insensitive  to  mesh  spacing.  The  pitot  pressure 
is  somewhat  sensitive  to  grid  refinement  as  the  streamwise  Mach  number  in  the  vortex  core 
increases  with  finer  meshes.  This  is  consistent  with  the  computed  vortex  core  sise  becoming 
smaller  with  mesh  refinement.  In  general,  the  solution  on  the  128  x  128  grid  can  be  considered 
converged  except  for  the  sise  and  details  of  the  inner  vortex  core  region[13].  However,  its  position 
and  maximum  total  pressure  loss  are  converged  as  illustrated  in  the  table.  The  conclusions  are 
representative  of  all  our  conical  Euler  calculations  for  supersonic  freestream  flows. 

The  second  example  of  a  conical  Euler  solution  is  taken  from  [15]  which  presents  six  different 
comparisons  between  conical  Euler,  conical  Navier-Stokes,  and  experimental  results.  Details  of 
the  calculations  and  comparisons  may  be  found  in  that  paper.  The  particular  example  presented 
here  is  for  a  delta  wing  with  a  leading  edge  sweep  A  =  75°  in  a  M»  =  2.8  flow  at  an  angle  of 


4P«rsonal  communication  between  authors. 


attack  a  =  8°.  These  parameters  give  Mn  =  .82,  a„  =  28.5°  which  lie  in  the  left  middle  region 
of  Figure  2. 
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Figure  9:  Total  pressure  loss,  Euler.  Figure  10:  Total  pressure  loss,  Navier-Stokes. 


Figure  11:  Cp,  Euler  and  data.  Figure  12:  Cp,  Navier-Stokes  and  data. 

A  comparison  of  the  total  pressure  loss  contours  between  the  conical  Euler  and  conical 
Navier-Stokes  (laminar)  is  given  by  Figures  9  and  10.  The  Navier-Stokes  model  clearly  shows 
the  viscous  boundary  layer,  the  secondary  separation,  and  for  this  case  a  tertiary  separation. 
The  latter  affect  the  position  and  shape  of  the  primary  vortex  by  their  displacement  effects. 
However,  both  the  Euler  and  Navier-Stokes  models  yield  a  primary  vortex  with  essentially  the 
same  total  pressure  loss  in  the  core.  A  comparison  of  the  surface  pressures  given  in  Figures  11 
and  12  shows  that  both  models  are  in  reasonable  agreement  with  the  data.  The  effects  of  the 
secondary  and  tertiary  vortex  (including  their  displacement  effects  on  the  primary  vortex)  are 
small,  but  well  represented  by  the  Navier-Stokes  solution. 


Total  Pressure  Loss 

The  origin  of  the  total  pressure  loss  in  the  Euler  equation  solutions  is  not  immediately 
apparent  and  has  been  the  source  of  considerable  controversy.  The  total  pressure  loss  can 


equivalently  be  represented  as  entropy  (5)  production,  and  by  Crocco’s  equation 

TVS  +  5xu  =  V/i,  (13) 

can  be  related  to  a  vorticity  (w)  distribution  since  the  total  enthalpy  (h0)  is  constant  for  this 
class  of  flows.  One  source  of  total  pressure  loss  is  the  conical  bow  shock  wave.  However,  these 
losses  are  two  to  three  orders  of  magnitude  smaller  than  the  values  in  the  vortical  regions. 
The  results  shown  in  Figures  3  to  8  contain  strong  embedded  shocks  which  clearly  have  an 
associated  total  pressure  loss.  Comparing  Figures  3  and  4  it  can  be  seen  that  these  losses 
are  approximately  convected  along  streamlines.  There  are  additional  total  pressure  losses  in 
the  shear  layer  from  the  wing  tip  and  vortex  core.  The  latter  is  clearly  not  convected  along 
streamlines.  Many  examples  have  been  computed  which  have  either  no  embedded  shocks  or 
very  weak  ones,  but  which  have  total  pressure  losses  of  30%  or  more.  It  is  clear  that  some 
mechanism  other  than  shock  waves  is  involved. 

A  possible  explanation  could  be  that  the  losses  represent  numerical  error  introduced  by  the 
discretization  process  or  added  damping  terms  as  represented  by  Equation  (12).  In  fact  this 
seems  to  be  true,  but  the  interpretation  is  a  bit  surprising.  Careful  numerical  studies  have  been 
done  [13,21]  in  which  grid  density  and  artificial  viscosity  constants  vj  and  were  varied  by  an 
order  of  magnitude.  The  total  pressure  loss  is  extremely  insensitive  to  these  parameters.  See, 
for  example,  the  table  in  the  previous  section.  The  total  pressure  loes  was  also  found  to  be 
insensitive  to  grid  distribution,  initial  conditions,  and  convergence  tolerance. 


Figure  13:  Second  order  damping  magnitude  Figure  14:  Equivalent  Reynolds  number 

A  clue  to  the  source  of  the  losses  can  be  found  by  comparing  the  magnitude  of  the  second 
order  damping  terms  with  the  sink  term  H  in  Equation  12  as  shown  in  Figure  13  taken  from 
[18].  The  results  are  for  a  A  =  75°  flat  plate  delta  wing  at  Af*  =  1.3,  a  =  10°  with  a  128  x  128 
cell  grid.  It  can  be  seen  that  in  the  feeding  sheet  region,  the  vortex  core,  and  around  an 
embedded  shock  under  the  vortex,  the  damping  terms  are  significant. 

Further  insight  can  be  obtained  by  looking  at  a  simplified  equation  set  corresponding  to  the 
one-dimensional  continuity  and  momentum  equations  [13] 

d  ,  ,  a* 

£<'*>  =  1,350  (U) 


(15) 


where  the  numerical  diffusion  coefficients  have  been  denoted  by  e.  Differentiating  momentum 
and  using  continuity  yields 


du  dp  .  d2u  _  *  d3p  dpd u 

pUdi+di~pefud^  +  u{€fU  tf)  dx'+^TxTz 


(16) 


Now,  if  e^t  =  if,  i.e.  the  same  artificial  viscosity  formulation  is  used  for  both  continuity  and 
momentum  equations  which  is  the  usual  procedure,  then  there  are  only  two  added  terms.  The 
first  is  a  diffusion  term,  while  the  last  one  is  a  convective  term  modifying  the  coefficient  of 
Powell  [13]  used  this  idea  to  define  an  effective  numerical  Reynolds  number  ~  i^}  which  varies 
throughout  the  flow  field  due  to  mesh  resolution,  the  magnitude  of  the  pressure  switch,  and  the 
value  of  A.  Figure  14  taken  from  [13]  shows  the  magnitude  of  the  numerical  Reynolds  number 
for  the  same  wing  as  Figure  13  but  at  M*  =  1.1,  a  =  10°.  An  embedded  grid  was  used  such 
that  in  the  region  of  interest  the  grid  density  is  256  x  256.  Even  for  such  a  fine  grid,  it  can  be 
seen  that  the  effective  Reynolds  number  in  the  primary  vortex  core  is  about  10s! 

It  is  worth  keeping  in  mind  that  most  Navier-Stokes  algorithms  contain  artificial  diffusion 
effects  (whether  added  explicitly  as  herein  or  included  implicitly  as  with  upwind  methods) 
which  are  needed  for  the  Euler  part  of  the  algorithm.  Unless  the  grid  is  extremely  fine  in  the 
vortex  region,  the  artificial  diffusion  terms  probably  dominate  the  real  viscous  terms  in  such 
calculations. 

Another  interesting  numerical  result  can  be  understood  from  Equation  16.  A  calculation 
was  done  in  which  the  second  order  damping  was  suppressed  from  the  continuity  equation 
(if  =  0)  but  retained  in  all  the  other  equations.  This  had  a  strong  influence  on  the  maximum 
total  pressure  loss.  It  can  be  seen  from  (16)  that  the  coefficient  i^  -  if  does  not  zero  out  and  an 
added  diffusion  term  results.  By  explicitly  correcting  for  this,  the  calculations  with  suppressed 
second  differencing  damping  in  the  continuity  equation  gave  the  same  total  pressure  loss  as  the 
original  computations.  So  the  magnitude  of  the  total  pressure  loss  is  sensitive  to  the  form  of 
the  added  viscosity,  but  not  to  its  magnitude.  Comparisons  between  the  present  finite  volume 
algorithm  and  results  from  other  algorithms  have  not  shown  significant  differences,  except  for 
the  results  of  Marconi  [22]  discussed  in  the  last  paragraph  of  this  section. 

Computational  results  indicate  that  the  maximum  total  pressure  loss  is  directly  related  to 
the  aerodynamic  parameters.  Variations  of  Mach  number,  angles  of  attack,  sweep,  and  yaw,  as 
well  as  body  geometry  all  affect  the  total  computed  pressure  loss  in  a  consistent  way.  That  is, 
increasing  the  angle  of  attack  produces  stronger  vortices  and  greater  total  pressure  loss  as  is 
illustrated  in  the  next  section. 


The  above  evidence  can  be  summarized  as  follows.  In  the  primary  vortex  core  and  the 
shear  layer  emanating  from  the  wing  tip,  the  added  artificial  viscosity  terms  are  significant 
and  produce  a  total  pressure  loss  in  a  manner  similar  to  what  the  actual  fluid  viscosity  would 
produce.  The  intriguing  question  is  why  the  total  pressure  loss  is  insensitive  to  the  magnitude  of 
the  numerical  diffusion  coefficient.  There  are  different  explanations  for  the  core  and  shear  layers. 
For  the  core  region,  the  circumferential  velocity  is  brought  to  zero  (see  Figure  7)  producing 
a  loss.  Powell’s  [13,23]  analytical  model  for  a  high  Reynolds  number  conical  incompressible 
viscous  vortex  shows  that  the  only  parameter  which  sets  the  maximum  total  pressure  loss  is  the 
swirl  velocity  at  the  edge  of  the  core.  As  can  be  seen  from  the  table  in  the  previous  section,  this 
parameter  is  also  independent  of  mesh  spacing  for  the  supersonic  conical  flows.  For  the  shear 
layer  [21],  the  added  viscous  terms  together  with  the  discrete  nature  of  the  solution  produce  a 
vortex  "sheet”  with  a  thickness  related  to  the  mesh  spacing.  The  velocity  magnitude  should 
be  the  same  on  either  side  of  the  sheet,  but  the  direction  changes  as  dictated  by  the  strength 
of  the  sheet.  TVaversing  across  the  sheet,  the  velocity  magnitude  will  have  a  minimum  in  the 
center  leading  to  a  total  pressure  loss.  The  reader  is  referred  to  [13,23,21,24]  for  the  detailed 
discussion  of  these  arguments. 


The  question  remains  as  to  whether  the  artificially  introduced  total  pressure  loss  is  related 
to  the  real  total  pressure  loss.  Comparisons  with  measured  pitot  pressure  as  given  in  the 
previous  section  and  total  pressure  as  given  in  the  next  section  indicate  that  although  the 
numerical  losses  may  not  be  exact,  they  are  a  good  approximation  for  compressible  flows.  For 
incompressible  flows,  the  agreement  is  not  as  good  [25],  but  the  studies  are  not  as  thorough. 
The  artificial  diffusion  terms  are  similar  to,  but  different  than,  the  real  diffusion  terms  even 
for  laminar  flows.  For  Reynolds  numbers  characteristic  of  aeronautical  applications,  the  vortex 
cores  will  be  turbulent.  A  more  physically  based  numerical  viscosity  model  would  be  helpful  in 
establishing  the  accuracy  of  the  predicted  total  pressure  losses. 

The  final  remark  of  this  section  addresses  the  issue  of  isentropic  vortex  cores.  Marconi 
[22]  has  suggested  that  the  correct  solution  of  the  Euler  equations  should  only  have  shock 
induced  total  pressure  losses.  He  has  used  the  A-scheme  which  integrates  the  equations  in  a 
charactersitic  form.  Entropy  production  must  be  introduced  explicitly  which  only  occurs  at 
shock  waves.  Vortical  flows  have  been  computed  with  zero  total  pressure  losses.  By  Crocco’s 
equation  (13)  this  means  the  vorticity  and  velocity  vectors  are  aligned.  The  results  of  these 
calculations  show  that  the  vortex  position,  surface  pressures,  and  other  global  features  of  the 
flow  are  in  close  agreement  with  solutions  which  have  total  pressure  losses  [22].  Powell  [13] 
modified  the  present  solution  algorithm  by  overriding  the  streamwise  momentum  equation  with 
the  constant  total  pressure  requirement.  Although  the  iterative  convergence  rate  is  poorer, 
solutions  can  be  obtained  and  compared  with  the  “lossy”  ones.  In  agreement  with  Marconi,  the 
global  features  are  the  same,  but  the  values  of  parameters  in  the  vortex  cores  are  different.  In 
particular,  the  streamwise  velocity  and  Mach  number  are  greater.  Comparison  with  the  pitot 
pressure  measurements  of  Monnerie  and  Werll  [16]  shows  much  poorer  agreement.  The  present 
authors  feel  that  the  desirable  solution  of  the  Euler  equations  is  that  which  models  a  viscous 
flow  in  the  limit  of  Re  -*  oo.  The  solutions  with  total  pressure  losses  due  to  artificial  viscosity 
are  more  realistic  than  those  with  no  total  pressure  loss.  Either  one  seems  to  be  acceptable 
for  prediction  of  surface  pressures  and  hence  forces.  However,  the  two  approaches  may  give 
different  results  regarding  the  prediction  of  vortex  bursting  in  three-dimensional  calculations 
as  total  pressure  losses  can  play  a  role  in  this  phenomenon  [26] . 


Three-Dimensional  Solutions 

Although  the  conical  Euler  model  is  suitable  for  studying  some  of  the  basic  features  of 
leading  edge  vortex  flows,  a  fully  three-dimensional  model  is  needed  for  realistic  wing  geometries 
as  well  as  for  subsonic  fireestream  Mach  number  conditions.  In  this  section,  some  results  using 
two  different  grid  topologies  are  presented  and  compared  to  experimental  data  for  Af,*,  —  .85. 
The  grids  are  rather  coarse  and,  therefore,  the  results  should  be  considered  as  exploratory.  The 
calculations  are  for  the  A  =  65s  sharp  edge  cropped  delta  wing  geometry  discussed  by  Elsenaar 
and  Biitefisch  [7],  Much  more  refined  grid  solutions  for  this  wing  are  presented  by  Wagner  and 
Hitzel  [27]  and  Kumar  and  Das  [28]. 

The  first  set  of  calculations  were  done  by  Malecki  [29]  using  an  0-0  grid  supplied  by  Rizzi 
and  his  co-workers  (e.g.  [8]).  A  total  of  37,768  computational  cells  were  in  the  field  with  2048 
cells  adjacent  to  the  wing  surface.  Although  1/16  of  the  total  cells  are  on  the  wing,  there  are 
only  16  cells  between  wing  and  far  field.  From  the  conical  Euler  studies,  this  is  not  adequate  to 
give  good  resolution  of  the  vortex.  Also  the  0-0  grid  has  significant  skewness  near  the  cropped 
tip.  Convergence  difficulties  from  this  were  solved  by  raising  v*  to  .075  in  this  area.  Only  10° 
angle  of  attack  was  reported  by  Malecki. 

The  second  set  of  calculations  were  done  by  Goodsell  [25]  who  modified  Roberts’  program 
to  handle  an  O-H  grid.  In  planes  perpendicular  to  the  flow,  an  O  grid  is  generated  similar  to 


th&t  used  for  the  conical  Euler  calculations.  A  total  of  73,728  cells  with  768  oq  the  wing  surface 
were  used.  There  are  only  12  chordwise  stations  on  the  wing,  but  there  are  48  cells  between  the 
wing  and  the  far  field.  Since  the  vortical  flow  does  not  vary  rapidly  in  the  streamwise  direction, 
this  grid  provides  better  resolution  of  the  crossflow  effects  than  the  0-0  grid.  However,  the 
streamwise  transonic  effects  such  as  shocks  will  be  poorly  resolved.  The  0-H  grid  has  less  skew¬ 
ness  than  the  0-0  grid,  and  has  a  topology  which  is  more  flexible  for  complicated  geometries. 
Goodsell  reported  results  for  a  =  10°,  20°, 24°,  and  the  former  two  are  discussed  here. 


Figure  15:  Surface  Mach  Number,  0-0  grid.  Figure  16:  Surface  Mach  Number,  0-H  grid. 
a  =  10°,Afoo  =  .85.  a  =  10°,  Moo  =  .85. 

The  results  for  a  =  10°  correspond  to  M«  =  .38,  on  =  22.6°  placing  them  well  within  the  left 
hand  middle  region  of  Figure  2.  Leeside  surface  Mach  number  contours  for  the  two  topologies 
are  shown  in  Figures  15  and  16.  The  footprint  of  the  vortex  is  evident,  and  the  flow  underneath 
it  is  supersonic.  The  0-0  grid  gives  a  maximum  surface  Mach  number  of  1.32  while  the  0-H 
grid  with  better  resolution  off  the  wing  surface  reaches  1.60.  The  flow  is  just  supersonic  over 
the  mid  portion  of  the  wing  demonstrating  slight  streamwise  transonic  effects  for  this  lowest 
angle  of  attack.  Finally,  both  solutions  show  a  shock  wave  at  the  tapered  tip  which  is  aligned 
with  the  streamwise  direction.  The  supersonic  outboard  flow  under  the  vortex  passes  through 
this  shock  to  turn  downstream. 

Figures  17-19  compare  the  computed  surface  pressure  coefficients  for  the  0-H  grid  (lines) 
with  the  measurements  (symbols)  at  30%,  60%,  and  80%  chord  locations.  The  windward  side 
agreement  is  good  except  at  30%  chord  where  the  centerbody  may  have  an  effect.  On  the 
leeward  side,  it  is  seen  that  the  calculations  predict  significantly  more  suction  and  that  the 
vortex  is  located  further  outboard  compared  to  the  experiments.  In  general  the  agreement  is 
not  too  good.  Figure  20  shows  similar  results  for  the  0-0  grid  at  80%  chord  location.  These 
calculations  show  leas  of  a  suction  peak  than  the  0-H  grid  results  and  in  general  have  poorer 
resolution.  The  0-0  grid  calculations  have  been  compared  with  Rizzi’s  code  using  exactly  the 
same  grid  and  are  in  good  agreement  [29]. 

Total  pressure  loss  contours  for  the  0-H  grid  calculations  at  78%  chord  are  shown  in  Figure 
21.  Experimental  [30]  total  pressure  contours  at  the  60%  chord  location  are  shown  in  Figure 
22  for  a  wing  with  the  same  planform  but  fiat  upper  surface.  The  computed  results  show  a 
loss  of  38%  while  the  measurements  indicate  a  loss  of  30-40%  which  is  in  qualitative  agreement. 
Malecki’s  calculations  with  the  0-0  grid  predicted  44%  losses  at  about  the  90%  chord  location. 

The  results  for  a  =  20®  correspond  to  Mn  =  .45,  =  40.7°,  again  in  the  left  hand  middle 

region  of  Figure  2.  All  results  are  from  the  0-H  grid  study  [25],  and  a  value  of  yj  =  .1  was 


Figure  23:  Surface  Mach  number,  O-H  grid.  Figure'24:  Surface  Cp  at  60%  chord,  O-H  grid, 
a  =  20°,  Moo  =  85.  a  =  20°,  =  .85. 


needed  for  convergence.  The  surface  Mach  number  contours  of  Figure  23  show  that  the  flow  is 
highly  transonic.  A  Mach  number  of  2.4  is  reached  under  the  vortex,  and  the  streamwise  flow 
in  the  center  of  the  wing  is  mostly  supersonic.  The  more  refined  calculations  of  Wagner  and 
Hitzel  [27]  show  that  this  supersonic  region  is  terminated  by  a  cross  flow  shock  which  interacts 
with  the  vortex.  The  present  grid  is  too  coarse  to  resolve  this.  The  clustering  of  the  contours 
indicates  that  there  is  a  shock  wave  under  the  vortex.  The  surface  Cp  results  at  the  60%  chord, 
Figure  24  location  likewise  indicate  a  strong  suction  peak  terminated  by  a  shock  wave.  The 
agreement  with  the  measured  Cp  data  is  poor  on  the  leeside  due  to  the  viscous  effects.  The 
secondary  vortex  and  shock  wave-boundary  layer  interactions  significantly  reduce  the  suction 
peak.  A  Navier-Stokes  model  is  needed  for  accurate  surface  pressure  predictions  in  this  speed 
regime.  The  computed  and  measured  total  pressure  contours  are  shown  in  Figures  25  and  26 
respectively.  Again,  there  is  reasonable  agreement  between  the  computed  loss  of  56%  and  the 
measured  value  of  ~  60%. 
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Figure  25:  (1  -  Pt/Pt »)  at  78%  chord,  O-H  Figure  26:  (Pt/Ptoo)  at  60%  chord,  measured, 
grid,  a  =  20°,  M«  =  .85.  a  =  20.5°,  =  .85. 


Conclusions 


Leading  edge  vortices  introduce  transonic  crossflow  effects  due  to  their  accompanying  high 
swirl  velocities.  Based  upon  both  calculations  and  experiments,  embedded  transonic  zones  and 
shock  waves  are  the  norm  rather  than  the  exception.  These  effects  are  in  addition  to  the  classical 
transonic  Sow  features  of  lifting  wings.  Numerical  solutions  of  the  Euler  equations  for  sharp  edge 
wings  appear  to  model  the  main  effects  of  the  primary  vortex.  Both  analytical  and  numerical 
evidence  point  to  the  artificial  viscosity  terms  as  the  source  of  vorticity.  The  numerical  damping 
models  routinely  used  by  most  algorithms  seem  to  lead  to  realistic  total  pressure  losses.  More 
physically  based  artificial  viscosity  models  would  be  desirable.  Total  pressure  deficits  may 
be  important  in  the  prediction  of  vortex  bursting.  For  transonic  freestream  Mach  numbers, 
the  influences  of  secondary  vortex  formation  and  shock  wave-boundary  layer  interaction  are 
important  and  outside  the  capabilities  of  an  Euler  equation  model.  For  supersonic  freestream 
conditions,  the  secondary  vortex  effects  have  less  influence  on  the  leeside  pressure  distribution. 
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